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Abstract. We reexamined the experimental evidences for the possible existence of 
the superconducting (SC) gap nodes in the three most suspected Fe-pnictide SC 
compounds: LaFePO, BaFe2(Aso.67Po.33)2, and KFe2As2. We showed that while the T- 
linear temperature dependence of the penetration depth A(r) of these three compounds 
indicate extremely clean nodal gap superconductors, the thermal conductivity data 
\iuYT,H^o i'^s[H,T)/T unambiguously showed that LaFePO and BaFe2(Aso.67Po.33)2 
are extremely dirty, while KFe2As2 can possibly be clean. This apparently conflicting 
experimental data casts a serious doubt on the nodal gap possibility on LaFePO and 
BaFe2(Aso.67Po.33)2- 
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1. Introduction 

Despite the intensive research effort since the discovery of the Fe-based superconductors, 
[Ij the pairing symmetry of this new class of superconducting (SC) compounds has 
not yet been settled. Early theories and experiments appear to best support the sign 
changing s-wave pairing state (denoted as s±-, s^-, or ±s-state in the literature). 
[21 m m ini |6] However, there exist several Fe-pnictide compounds that are not seemingly 
compatible with the ±s state but are strongly suggesting for the presence of nodes in 
their SC states. Among others, LaFePO, [71 |8] BaFe2(Aso.67Po.33)25[9] and KFe2As2 
[To] are the most compelling compounds for the nodal gap, and in lesser degree 
Ba(Fei_2;Co2:)2As2 [llj is also suspected. 

Commonly taken evidences for the nodal gap in the above mentioned compounds 
are: (1) T-linear temperature dependence of penetration depth A(T) down to very low 
temperatures. [71 l8l l9l [TO] and (2) a strong field dependence in the thermal conductivity 
slope liuiT^o K,{H,T)/T, which is proportional to in between \/H and H, accompanied 
with a substantial fraction of the residual thermal conductivity limT^n^o ^^iHjT) /T = 
i^so/T. [01 [121 nil [13] These features are the well known signatures of the nodal 
gap superconductors such as the rf-wave superconductivity of the high-Tc cuprates. 
And although it was recently shown that the strong field dependence of the thermal 
conductivity k{H,T — )■ 0) can be equally well explained with the ±s-wave state, [U] 
the extremely close T-linear A(T) is hard to be reconciled with other than a clean nodal 
gap superconductor. Furthermore, the finite value of the residual thermal conductivity 
Kso/T measured in all three compounds [9l HH [131 E] ^it is known that the nodal gap 
SC state produces an universal thermal conductivity slope independent of the amount 
of impurity concentrations [161 113 IlH] - is another evidence for a nodal gap state, so it 
was widely interpreted to support the presence of nodes in these compounds together 
with the penetration depth data. 

In this paper, however, we will show that there is a serious and unreconcilable 
conflict between the above mentioned two experimental evidences for the nodal gap. 
We notice that (1) the universal value of ksq/T delivers no information about the 
dirtiness of the superconducting sample, however, (2) the normal state value of Kn/T 
tells us the amount of dirts in the sample. Then combining the facts (1) and (2), the 
ratio Ks/i^N, which is usually plotted data in experiments, is a very good indicator of 
the dirtiness of the sample. Inspecting the reported data of thermal conductivities 
of the three compounds, we concluded that the measured samples of LaFePO and 
BaFe2(Aso.67Po.33)2 should have a large amount of impurities and hence cannot be 
compatible of the T-linear A(T) within the nodal gap scenario. In the case of KFe2As2, 
there exist two very different thermal conductivity data by Dong et al. [13] and Reid et 
a/.[T^ with different Tc samples, 3K and 3.8K, respectively. Our analysis of the thermal 
conductivity data showed that the sample of [14J (T^ ~ 3.8K) is cleaner with at least 
10 times less impurity concentration than the sample of Ref. [T3] (Tc ~ 3K). Hence, 
the former sample of KFe2As2 can possibly be compatible with the T-linear A(T) data. 
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We conclude that KFe2As2 can remain a possible nodal gap superconductor, but not 
LaFePO and BaFe2(Aso.67Po.33)2- 



2. Theory 

Assuming the quasiparticle excitation E(k) = a/ (fi/ci)^ + (^^^2)^ in the d-wave 
superconductor {vp Fermi velocity perpendicular to the Fermi surface(FS), Vi nodal 
gap velocity parallel to the FS), the universal thermal (electric) conductivity in the 
nodal gap superconductor has been derived as follows, [ini [13 HB] 



(1) 







where Aq is the maximum gap value of the (i-wave gap A(0) and 7^ is the impurity 
induced damping rate at zero energy in the SC state. As well known, ^ indeed becomes 
universal, independent of the impurity concentrations and scattering strength, but only 
in the limit of Aq 3> 7^; Eq.(l) clearly shows that a deviation occurs when 7^ ~ Aq. 
The normal state limit of the above is easy to derive as 



/£iv(T = 0) vl Ao 

' \^) 

where 7„ is the impurity induced damping rate in the normal state and in general 7^ 7^ 7^ 
for the same impurity strength and concentration. Knowing that the normal state should 
have no memory of the superconductivity, the above expression of ^ is a disguised form 
for convenient comparison with Eq.(l) and ^ becomes 7r/4 or a momentum scale of 
the FS size. Also we don't need to know the material specific parameters like iV(0),fi?, 
etc to estimate the absolute magnitude of the thermal conductivities because for our 
purpose we only need the ratio 



r 1 r 

= . ^ — . (4) 

Ao v^l + r/Ao Ao 

where F = nimp/T^N{Q) is the impurity concentration parameter and in the second line 
of the above equations we used the results of 7^ ~ -v/AoF and 7„ = F assuming the 
unitary impurity scattering strength. Eq.(4) is the key result of this paper. 

It was nice to observe the universal value of the thermal conductivity slope k^/T of 
Eq.(l) to confirm a nodal superconductor. On the other hand, it was also a drawback 
since the universal thermal conductivity slope doesn't tell us how dirty or clean the 
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sample is. However, as shown in Eq.(4), the ratio ks/i^n is an excellent indicator of the 
dirtiness of the specific SC sample. At this point, we would like to recall the fact that 
the typical experimental values of ns/f^N measured in LaFePO and BaFe2(Aso.67Po.33)2 
are ~ 0.2 — 0.4 which is quite high level of impurity concentration for a nodal gap 
superconductor. 



3. Numerical calculations and discussions 



In this section we will show the full numerical calculations of the field dependence of 
lim-r^o i^siH, T)/T as well as the specific coefficient limx-^Q C{H, T)/T of the canonical 
d-wave gap (A(0) = Aocos2^) state with the various impurity concentrations. Also 
the results of the penetration depth A(T) will be shown for the corresponding impurity 
concentrations. 

To calculate the field dependencies of the thermal conductivity and specific heat 
in the mixed state with applied field, we just need to calculate the position dependent 
DOS A^(ci;, r) in the presence of vortices. In the semiclassical approximation, the matrix 
form of the single-particle Green's function in the SC state, including Doppler shift of 
the quasiparticle excitations e{k) due to the circulating supercurrent Vs(r), is given by 
[T511I91120] 

G(k, r, 6) = ^"'"-^t.''^I;'"'^tM^Ar (5) 
[w + Vs(r)-k]2-e2(A;)- A2 ^ ' 

where Tj are Pauli matrices and the supercurrent velocity around the vortex core 
Vs(r) is ~ mr ' "^i^h r the distance from the vortex core. The position dependent 
DOS is calculated as N{oj,r) = — ;iTrIm^^ J d9 Go{k,r,uj,9). Finally, the field 
dependent quantities are obtained from the areal average DOS per unit volume as 
N{uj,H) = J^" dr^N{uj,r)/7rR'jj with the magnetic length Rh = ($o a fiux 

quanta) and the SC coherence length ^. 

The impurity scattering is included by the T-matrix method. [2T1 123] The 
impurity induced self-energies renormalize the frequency and order parameter (OP) as 
uj ^ u = uj + TP{uj) and Aq -> Aq = Aq + ^^(w), with SO'i(w) = F ■ T^'H^^), where r°'^ 
are the Pauli matrices r°'^ components of the T-matrices in the Nambu space. However, 
is identically zero in the ci-wave state. Then all impurity effect and the Volovik effect 
can be incorporated into the local Green's function Eq.(5) by replacing uj by Co. 

After calculating the averaged iV(a;, H) for all frequencies, specific heat is calculated 

as 

CiT,H)=rM^f^^^^ (6) 
la t' cosh"(g.) ^ ' 

Similarly, thermal conductivity is calculated by [21] 

n(T,H,r)<.v,JdgJ^ (7) 
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Figure 1. (Color online) Thermal conductivity ks(H)/T vs the normalized fields 
H/Hc2 of the c?-wave SC state, calculated at T = 0.02Ao for various impurity 
concentrations F/Aq = 0.01, 0.02, 0.05, 0.1, 0.2, and 0.4. (unitary impurity). The inset 
shows the full range of fields up to H/ Hc2 = 1- 
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where z = a; + Vs(r) • kp. And the longitudinal and transversal thermal 
conductivities are calculated cLS K\\ 

j^" d^rK-\T, H, r)/7TRl, respectively. 
3.1. Thermal conductivity 

Figure 1 shows the theoretical thermal conductivity k,{T,H)/T vs H of the d-wave 
SC state calculated at the low temperature limit of T/Aq = 1/50 with the varying 
impurity concentrations of the unitary scatterers, F/Aq = 0.01,0.02,0.05,0.1,0.2, 
and 0.4. First, the results indeed showed that the universal thermal conductivity 
limT^H^o i^(T, H)/T is well reproduced by our numerical calculations for the vast 
range of impurity concentrations. Second, it showed that the normal state limit of 
K,{T, H)/T, which is approached by increasing the field strength H toward Hc2, is 
inversely proportional to the impurity concentration as shown in Eq.(4). The inset 
shows the results for the full range of H/Hc2 = [0:1] and we can see that n{H)/T 
sharply increases near Hc2- This is due to a rapid collapse of the gap Ao(-ff) towards 
Hc2 and our semiclassical approximation faithfully follows the Doppler shifting effect of 
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this rapidly collapsing gap up to Hc2- While this is the correct calculation results with 
the semiclassical approximation, it is also known that this semiclassical approximation is 
not precisely correct near Hc2 where the quantum effect should become important. [25] So 
the exact field dependence of k,{T, H)/T near Hc2 in Fig.l should not be taken seriously. 

However, the important points for our purpose are: (1) at both limits, the universal 
limit value of \im.T i^{T , H) /T and the normal state limit value n{H = Hc2)/T are 
exact, and (2) the overall field dependence of the initially slow rise and then a rapid rise 
of K,{H)/T near Hc2 is the genuinely correct behavior regardless of different theoretical 
treatments. [26] The main conclusions of this paper relies only on these two points. The 
main panel in Fig.l shows the results for the limited region of H < 0.9Hc2 for a better 
resolution of the low field behavior of k,s{H)/T. 

In Fig. 2, we replotted the theoretical results of Fig.l with two different 
normalizations and the experimental data of [ks{H)/T]/[k,n/T] are overlaid. Fig.2(a) 
normalized ks{H)/T of Fig.l with its ks{H = Hc2)/T values for each impurity 
concentration, and Fig. 2(b) used the ks{H = 0.9Hc2)/T values for normalization. The 
second normalization plot by the H = 0.9Hc2 values was chosen because it is the typical 
point of saturation before the sharp rise as seen in the inset of Fig.l. and this concave- 
down saturation behavior approaching Hc2 is the typical observations in experiments. 
[U [121 [13] Different normalizations yield different line shapes of the normalized ks{H)/T 
which is supposed to be compared to the experimental [Ks{H)/T]/[K]\r/T]. The true 
behavior should be somewhere in between Fig. 2 (a) and Fig.2(b), but we emphasize 
that this fine detail is irrelevant to our main conclusions and analysis. The overlayed 
experimental data are BaFe2(Aso.67Po.33)2, [Sj LaFePO, [12] and KFe2As2. [T3[ [T^ 

Regardless of the choice of the normalizations, the experimental values of the 
residual thermal conductivity of BaFe2(Aso.67Po.33)2, [S] LaFePO, [12] and KFe2As2 [13] 
unambiguously indicate that these compounds should have the impurity concentration 
F/Aq > 0.4, which is extremely dirty superconductor. Our theoretical calculations are 
with a single rf-wave gap band. In reality, if there exists a nodal gap in these multiband 
Fe-pnictide compounds, the total gap function should consist of a nodal gap + one or 
two s-wave gaps, for example, a nodal ±s-wave gap.[l2j If that is the case, the total 
ktv/T should increase due to the addition contributions from other bands. However, 
these additional s-wave gap bands have negligible contributions to the residual thermal 
conductivity limx^H^o f^s{H) /T because they are fully gapped at low fields and low 
temperatures. Therefore, we need to have even higher impurity concentration than 
F/Aq > 0.4 in order to match the experimental data[9l|12l[T3] of the normalized residual 
thermal conductivity \imH,T~^o[Ks{H) /T]/[kn/T]. 

On the other hand, the data of KFe2As2 by Reid et al. [13] is very different from 
the data of KFe2As2 by Dong et al. [13] We can see that the data of Reid et al. 
[H] reasonably fit the theoretical result in Fig. 2 (a) with the impurity concentration, 
F/Aq < 0.02, which is relatively clean limit. As discussed in Ref.[Tl], the discrepancy 
between the data of two groups is understood by the sample purity. Judging from the Tc 
of two samples (3.80K and 3K, respectively) and our theoretical calculations of hs/i^n 
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Figure 2. (Color online) (a) Theoretical thermal conductivities (lines) of Fig.l 
normalized by the values of ks{Hc2)/T for each impurity concentrations F/Aq = 
0.01, 0.02, 0.05, 0.1, 0.2, and 0.4. (the same colors as in Fig.l). Symbols are the cropped 
data of experiments, (b) The same data as in (a) but normalized by the values of 
Ks{H — 0.9Hc2)/T for each impurity concentration. Symbols are the same as in (a). 
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in Fig. 2, the sample of Reid et al. must be much cleaner, by about 10-20 times, than the 
one of Dong et al. and it appears to be consistent with the result of the clean d-wave 
calculation with F/Aq < 0.02 in Fig.2(a). 

Summarizing the cases of BaFe2(Aso.67Po.33)2 [9] and LaFePO [12], if we interpret 
the thermal conductivity data of these two compounds with a nodal gap scenario, we 
are led to conclude that both compounds are very dirty nodal gap superconductors. 
However, if this is true, the linear-T penetration depth measurements |7l El [9] of these 
two compounds are in a serious conflict with the dirty nodal gap scenario. 

In the case of KFe2As2, we have two seemingly contradicting thermal conductivity 
measurements [131 E] seen in Fig. 2 (a) and (b). However, despite the large difference 
of Tc and the line shapes of k,s{H)/ kn, both samples reported a similar value of the 
residual thermal conductivity: nso/T=3.7 ± 0.4 mW/K^cm (Ref. [T4]) and k,so/T=2.27 
± 0.02 mW/K^cm (Ref.[i3j, here we ignored the correction by geometric factor discussed 
in Ref. [H]), respectively. This fact itself is a strong supporting evidence for the nodal 
gap in KFe2As2 compound. The different zero field intercepts of the data of two samples 
in Fig. 2 are due to the normalization by the normal state thermal conductivities: 
K^/T=7.36 ± 0.04 mW/K^cm for Dong et al. [13] and kn/T ^ 109 mW/K^cm for 
Reid et al. [H] The data of Reid et al. [14J can fit reasonably well with the clean nodal 
gap calculation (F/Aq < 0.02) for the most of low field region as seen in Fig.2(a). The 
data of Dong et a/.,[l3] however, doesn't fit with any calculational results in Fig.2. 
While the residual thermal conductivity value of it can be fit with a dirty nodal gap 
with F/Aq ~ 0.4, the data for H > increases much rapidly at low fields and saturates 
to become fiat for H > 0.4i7c2; not even close to any theoretical results in Fig.2. But, 
if we assume the additional bands with s-wave gaps in addition to a nodal gap, as in a 
nodal ±s-wave state, it would be possible to fit the data of both clean and dirty limits 
of Ref. [131 [H]- The detailed SC properties of the nodal ±s-wave state will be reported 
in future publication. 

3.2. Specific heat coefficient and Superfluid density 

To foster the above discussions, we calculated the field dependence of the specific 
heat coefficient '~is{,H) = \imT^QC{H,T)/T and the temperature dependence of the 
superfluid density p(T) ~ 1/A^(T) of the d-wave state with various concentrations of 
the unitary impurities as in Figs.l and 2. Figure 3 is the normalized 7<j(i?)/7n for 
F/Aq = 0.4, 0.2, 0.1, 0.05, 0.02, and (no impurity limit). It shows the expected behavior 
as the y/H behavior of 7s (-ff) in the clean limit becomes flattened with increasing the 
impurity concentration. The only point that we want to emphasize for our purpose 
is that even a small amount of impurities, for example, F/Aq = 0.02, immediately 
creates a substantial fraction of the specific heat coefficient 7s (if = 0) ~ 0.27^. This 
demonstrates that a nodal gap such as the d-wave state is extremely vulnerable to the 
unitary impurity scattering to create the low energy excitations. 

Figure 4 shows the normalized superfluid density ps{T) ~ 1/A^(T) of the rf-wave 
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Figure 3. (Color online) Normalized specific heat coefRcient lim-T^o C{T,H)/T = 
75 (i?) vs fields H/Hc2 calculated at T = 0.02Ao of the d-wave SC state for various 
impurity concentrations F/Aq = 0,0.02,0.05,0.1,0.2, and 0.4 (unitary impurity). 




Figure 4. (Color online) Normalized superfluid density ps{T) vs T/Tc of the d-wave 
SC state for various impurity concentrations F/Ao — 0,0.02,0.05,0.1,0.2, and 0.4 
(unitary impurity). 2Aq/Tc = 4 is used. 
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state with corresponding impurity concentrations of Fig.3. It also shows the well 
known behavior of ps{T) of the d-wave state with impurities. The typical T-linear 
ps{T) behavior in the clean limit changes to the behavior at low temperatures with 
impurities. Similarly to the evolution of the specific heat coefficients in Fig.3, even 
a small amount of impurities changes quite a wide temperature region into the 
behavior. For example, the impurity concentration of F/Aq = 0.02 makes ps{T) ~ 
for < T < 0.2Tc. In view of the fact that all three pnictide compounds studied in this 
paper reported the T-linear behaviors of A(T) down to extremely low temperatures: 
LaFePO ( 0.02 < T/T^, [7] and 0.08 < T/T^ [8J), BaFe2(Aso.67Po.33)2 (0.025 < T/T^ 
[9]), and KFe2As2 (0.05 < T/Tc [lOj), the required purity of these samples for the 
nodal gap scenario is F/Aq < 0.01 or even cleaner. However, this clean nodal gap 
scenario is totally contradicting to the thermal conductivity measurements ks{H)/T in 
the cases of BaFe2(Aso.67Po.33)2 [2] and LaFePO, [12] but possibly not with the case of 
KFe2As2.[T3l[Tl] 

4. Remark on the quantum oscillation (QO) experiments 

It is well known that the observation of the QO such as de Haas-van Alphen (dHvA) 
oscillation is possible only with very clean samples and, therefore, often quoted as an 
indication of the extreme purity of the probed sample. All three compounds studied 
in this paper have reported the QO experiments [271 EHl EH], hence we should worry 
about the consistency between the QO experiments - which warrant that the probed 
samples are clean - and our analysis of the thermal conductivity - which indicate that 
most of theses compounds are not so clean. To begin with, we note that the criteria 
of the cleanness for the QO signal and the SC properties are different; the former is 
F < hwc{= ^C-B/m* ; m* =renormalized mass) and the latter is F < Aq. Also, we need 
some interpretation for the typical dirtiness deduced from our thermal conductivity 
analysis, i.e., F/Aq, which was based on the single (i-wave model. We concluded in 
the previous sections that only KFe2As2 is possibly consistent with a nodal gap, but 
LaFePO and BaFe2(Aso.67Po.33)2 are not consistent with a nodal gap state but would be 
more consistent with a ±s-wave gap model with a small isotropic gap on the major band 
and a larger isotropic gap on the minor band, namely, A5 ^ A^, and A^(0)5 ^ N{0)l. 
Therefore, in the cases of LaFePO and BaFe2(Aso.67Po.33)2, for example, the deduced 
damping rate F/Aq ~ 0.4 should be understood as F/A5 ~ 0.4 while T^. of the com- 
pound is mostly governed by Al.[30] Bearing these in mind, let us examine the cases of 
each compound below. 

4.1. LaFePO (T^ ~ QK). 

The QO measurement on LaFePO has observed signals with magnetic fields above ~ 9T 
but the practically useful signals were obtained above ~ 20T and up to 45T.[27j The 
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estimated cyclotron frequency at B = 20T is fkOc ~ 1-2 — 2AmeV with the renormahzed 
mass m*/mo (mo, free electron mass) ~ 1 — 2. [27] and the estimated damping rate 
from our analysis is F = 0.4Ao ~ 0.35mey assuming the BCS relation Aq/Tc = 1.75. 
Therefore, the condition for the QO observation fiWc > T sufficiently holds for all fields 
B > 20T, and, therefore, without invoking further argument of the multiple gaps. As 
and Al, the observation of the QO in LaFePO has no contradiction with the damping 
rate estimated in our analysis. 

4.2. KFe2As2 (T^ ~ 3K). 

The QO signals on KFe2As2 [28j were obtained in the field range of 10 to 17.5 T. 
The estimated cyclotron frequency is hjjc ~ 0.2meV aX B = lOT with the heavily 
renormalized mass m*/mo ~ 6 of this compound. |28j As discussed in the previous 
section, there exist two very different thermal conductivity experiments |1 3 [ ITi] : the 
sample by Reid et a/. |14j seems to be clean (F < 0.02Ao) and the other one by 
Dong et a/. [13] seems to be dirtier (F ^ 0.4Ao). The estimated damping rates are 
F ~ 0.0087me^ for the clean one and F ~ 0.175mel^ for the dirty one, respectively, 
assuming the BCS relation Aq/T^ = 1.75 with =3K. Therefore, if the sample used 
for the QO experiment [28] is close to the cleaner one, there is absolutely no problem 
to observe the QO signals {huc ^ F; hwc ~ 0.2meV, F ~ 0.0087meV^). On the other 
hand, if the sample were on the side of the dirtier one, the observation of the QO signals 
should be very weak at best (F ~ 0.175meV, hojc ~ 0.2meV). Putting together, there 
exists a wide range of sample purity between 0.0087meV" < F < 0.175mel^ with which 
the QO experiment was possible. 

4.3. BaFe2(Asos7Po.33)2 (T.c^^QK). 

Shishido et a/. [29] have performed the QO experiments with BaFe2(Asi_.a;P2:)2. However, 
the QO signals was obtained only with 1 > x > 0.41 for the field range from 17T to 
55T, and the x = 0.33 sample never produced meaningful signals up to 55T. On top of 
that, even in the samples of 1 > x > 0.41 only the electron band FSs {a and /3 bands 
in their notations) produced signals but the hole band FSs never produced measurable 
signals. With theses, we can estimate the overall damping rate of the x = 0.33 sample 
should be higher than hwc ~ 0.55meV using the renormalized mass ■m*/mo ~ 3 and 
the maximum field strength B = 55T used in experiments. [29] On the other hand, our 
estimated damping rate is F = 0.4Ao ~ l.imeV using the BCS relation and Tc=30K. So 
it is consistent with the failure of the QO experiment for the x = 0.33 sample. In reality, 
since we have argued that the ±s-wave state is more consistent with BaFe2(Aso.67Po.33)2, 
if we understood F = 0.4Ao as F = O.4A5, the real damping rate F should be < l.AmeV 
but still > 0.55me\^. 
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5. Conclusions 

In conclusion, we have carefully reexamined the experimental evidences for the possible 
existence of the SC gap nodes in the three most suspected Fe-pnictide compounds, 
LaFePO, BaFe2(Aso.67Po.33)2, and KFe2As2. We have derived an exact relation for the 
ratio between the universal residual thermal conductivity ks/T and its normal state 
value ktv/T in the d-wave state. Using this ratio ks/i^n ~ T/^o as an indicator to 
determine the dirtiness of the SC sample, we have shown that the reported experimental 
data of the thermal conductivity in BaFe2(Aso.67Po.33)2 [9] and LaFePO [I2] indicated 
that the measured samples are the dirty limit superconductors, hence contradicting to 
the clean limit nodal gap scenario deduced from the penetration depth measurements. 
[H [HI [9] To this end, if the nodal gap scenario fails for these two compounds, we propose 
a dirty ±s-wave gap state as a possible scenario to reconcile the apparently contradicting 
experiments of thermal conductivity and penetration depth measurements - the large 
residual thermal conductivity slope ks/T and the T-linear A(T); in this scenario one 
isotropic s-wave gap is much smaller than the other one and the small gap is almost 
filled with the impurity band caused by a sufficient amount of impurity scattering. 

In the case of KFe2As2, there exist two qualitatively different data of the thermal 
conductivity measurements. [131 113 appears that the one of Reid et al. is a clean 
sample but the one by Dong et a/. [13] contains at least 10 times more impurities. We 
concluded that the clean sample result can be consistent with a nodal gap scenario 
both for the thermal conductivity and penetration depth measurements. On the other 
hand, the thermal conductivity data of the dirty sample [13] can be understood with a 
dirty nodal gap plus additional isotropic gaps as in the nodal ±s-wave gap, but the T- 
linear A(T) cannot be compatible with this sample by any scenario. Therefore, the gap 
symmetry of KFe2As2 should be further investigated by the cross-examination of the 
penetration depth and thermal conductivity measurements with samples with various 
purities. 
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